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Abstract 

We propose that the Baxter’s Q-operator for the quantum XYZ spin chain with 
open boundary conditions is given by the j ^ oo limit of the corresponding transfer 
matrix with spin-j (i.e., {2j + l)-dimensional) auxiliary space. The associated T-Q 
relation is derived from the fusion hierarchy of the model. We use this relation to 
determine the Bethe Ansatz solution of the eigenvalues of the fundamental transfer 
matrix. The solution yields the complete spectrum of the Hamiltonian. 
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1 Introduction 


In this paper we are interested in the exact solution of the quantum XYZ quantum spin 
chain with most general nondiagonal boundary terms, dehned by the Hamiltonian 
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where a®, cr^ are the usual Pauli matrices, N is the number of spins, the bulk coupling 
constants J^, Jy, Jz are related to the crossing parameter rj and modulus parameter r by 
the following relations, 
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and the components of boundary magnetic helds associated with the left and right boundaries 
are given by 
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Here <7 (m) is the (7-function dehned in the next section and are free boundary param¬ 

eters which specify the boundary coupling (equivalently, the boundary magnetic helds). 

In solving the closed XYZ chain (with periodic boundary condition), Baxter constructed 
a Q-operator P], which has now been proved to be a fundamental object in the theory of 
exactly solvable models [2]. However, Baxter’s original construction of the Q-operator was 
ad hoc and its connection with the quantum inverse scattering method was not clear. It was 
later argued [S] that the Q-operator for the XXZ spin chain with periodic boundary condition 
may be constructed from the 14 ( 5 / 2 ) universal L-operator whose associated auxiliary space is 
taken as an inhnite dimensional q-oscillator representation of Uq{sl 2 )- (See also for recent 
discussions on the direct construction of the Q-operator of the closed chain.) In j3], a natural 
set-up to dehne the Baxter’s Q-operator within the quantum inverse scattering method was 
proposed for the open XXZ chain. However, the generalization of the construction to off- 
critical elliptic solvable models [H] (including the XYZ chain) is still lacking even for the 
closed chains. 
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In this paper, we generalize the construction in to Baxter’s eight-vertex model and 
propose that the Q-operator Q{u) for the XYZ quantum spin chain (closed or open) is given 
by the j —>■ cx) limit of the corresponding transfer matrix with spin-j (i.e., (2j -|- 1)- 

dimensional) auxiliary space, 

Q{u) = lim t^^\u — 2jiri). (1.3) 

jr'^OO 

This relation together with the fusion hierarchy of the transfer matrix leads to the T-Q 
relation. We then use the T-Q relation, together with some additional properties of the 
transfer matrix, to determine the complete set of the eigenvalues and Bethe Ansatz equations 
of the transfer matrix associated with the Hamiltonian (HID under certain constraint of the 
boundary parameters (see ()5.1f)|l below). 

The paper is organized as follows. In Section 2, we introduce our notation and some basic 
ingredients. In Section 3, we derive the T-Q relation from dn and the fusion hierarchy of 
the open XXZ chain. In section 4, some properties of the fundamental transfer matrix are 
obtained. By means of these properties and the T-Q relation, we in Section 5 determine the 
eigenvalues of the transfer matrix and the associated Bethe Ansatz equations, thus giving 
the complete spectrum of the Hamiltonian (HID under the constraint ()5.16|) of the boundary 
parameters. We summarize our conclusions in Section 6. Some detailed technical calculations 
are given in Appendices A-C. 


Fundamental transfer matrix 


Let us £x T such that Im(r) > 0 and a generic complex number p. Introduce the following 
elliptic functions 
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a{u) = 6 

Among them the cr-function^ satishes the following identity: 

a{u + x)a{u — x)a{y + y)a{y — y) — cr{u -I- y)a{u — y)a{y x)a{y — x) 

= a{u + v)a{u — v)a{xy)a{x — y), (2.3) 

^Our cr-function is the t?-function i?i(it) [ 7 |- It has the following relation with the Weierstrassian a- 
function if denoted it by aw{u): aw{u) oc <t{u), rji = 7r^(i — ^nd q = 


3 










which will be useful in deriving equations in the following. 

The well-known eight-vertex model R-matrix R{u) G End(C^ 0 C^) is given by 


R{u) = 

The non-vanishing matrix elements are [2j 
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Here u is the spectral parameter and rj is the so-called crossing parameter. The R-matrix 
satishes the quantum Yang-Baxter equation 


Ri,2(mi - U2)Ri,3{ui - U3)R2,3{u 2 - U3) = R2,3{u2 - M3)-Ri,3(mi “ M3)-Ri,2(mi “ U2), (2.6) 

and the properties, 

PT-symmetry : Ri, 2 {u) = R 2 ,i{u) = R^ 2 ^(-u), (2.7) 

Z2-symmetry : a\a''2Ri^2iu) = Ri^2iu)a\a\, ioTi = x,y,z, (2.8) 

Unitarity relation : 2 (m)R 2 i(— m) = —^(m) id, ^iu) = — —— 7 -^—(2.9) 

a{T])a{T]) 

Crossing relation : Ri, 2 (m) = ViR^ 2 (“'^ ~ vWi: ^ (2-10) 

Here i? 2 ,i(w) = Pi 2 Ri, 2 {u)Pi 2 with P 12 being the usual permutation operator and ti denotes 
transposition in the Tth space. Throughout this paper we adopt the standard notations: for 
any matrix A G End(C^), Aj is an embedding operator in the tensor space ® ® • • •, 

which acts as A on the j-th space and as identity on the other factor spaces; Ri,j{u) is an 
embedding operator of R-matrix in the tensor space, which acts as identity on the factor 
spaces except for the Tth and j-th ones. 


4 









































Integrable open spin chains can be constructed as follows [S]. Let us introduce a pair of 
K-matrices K~{u) and K^{u). The former satishes the reflection equation (RE) 


Ri, 2 (Mi - M 2 )R'i (Mi)R 2 ,i(Mi + M 2 )R '2 (“ 2 ) 

= {u2)Ri,2{ui + U2)K{{ui)R2,i{ui - M 2 ), ( 2 . 11 ) 

and the latter satishes the dual RE ElEI 

Ri,2{u 2 - Ui)K^{ui) R2 ,i{-Ui -U2- 27]) iu2) 

= K^{u2)Ri,2{-Ui -U2- 27])Kt{ui)R2,l{u2 - Ml). (2.12) 

Then the transfer matrix t{u) of the open XYZ chain with general integrable boundary terms 
is given by [S] 

t{u) = tro (^K+ {u)To{u)K^{u)fo{u)'^ , (2.13) 

where To(m) and To(m) are the monodromy matrices 

To{u) = Ro^n{u) ... i?o,i(w), To{u) = Rifi{u)... Rn,o{u), (2.14) 

and tro denotes trace over the “auxiliary space” 0 . 

In this paper, we consider the most general solutions K^{u) [TU] to the associated rehec- 
tion equation and its dual, 
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R+(m) 
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where / is the 2x2 identity matrix and the constants are expressed in terms of 

boundary parameters as follows: 
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Sklyanin has shown that the transfer matrices with different spectral parameters commnte 
with each other: [t{u ), t(n)] = 0. This ensnres the integrability of the open XYZ chain. The 
Hamiltonian (HH) can be expressed in terms of the transfer matrix 




(2.18) 


where a'(0) = 

3 Fusion hierarchy and T-Q relation 

We shall nse the fusion procedure, which was hrst developed for R-matrices ffH and then 
later generalized for K-matrices [HI ESI, to obtain the Baxter T-Q relation. The fused 
spin-(j, i) R-matrix (j = 1, |,...) is given by [TT] 



^(i...2i)(“) = Ptlj {K2jiu)R2j,2j-i{‘2u + r])K^j_^{u + r]) 


xR2j,2j-2i‘2u + 2r])R2j-i,2j-2{2u + 3ri)K2j_2{u + 2ri) ■ ■ ■ 

X R2,i{2u + (4j - 3)r])K{{u + {2j - I)?])} (3.2) 


The fused spin-j K-matrix given by 



(3.3) 


where the scalar functions F{u\q) are given by F{u\q) = 1/(HLi nl;=i 'C(2^ + Irj + kr])), for 


q = 1,2,.... The fused transfer matrix tR\u) constructed with a spin-j auxiliary space is 
given by 



where 


T{l,„2j){u) — P(1...2j),Ar(w) . . . P(1...2j),l(M), 
T{l...2j){u + (2j - 1)t]) = R{l...2j),l{u) . . . R{1...2j),N{u). 


(3.5) 
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The transfer matrix corresponds to the fnndamental case j = |, i-e., t{u) = 

The fnsed transfer matrices constitnte commntative families, namely, 

= 0. (3.6) 

They satisfy a so-called fnsion hierarchy relation [121 cni 

(3.7) 

In the above hierarchy, we have nsed the convention = id. The coefficient fnnction 

6{u) can be expressed in terms of the qnantnm determinants of the monodromy matrices 
T{u), T{u) and K-matrices |T21IT!^ . 


t^^\u-{2j - l)vi)=t'^^ 2 )(M-( 2 j - l)ri)t{u) - - l)r]), j = 1, -, 
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Here pi,i{u) = and the corresponding determinants are given by 
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Det{K“(u)} = tri2(^Pi2^K^{u-p)Ri2{2u-p)K2iu)^ , 

Det {K+(n)} = tru (^Pi 2 ^ K^{u)Ru{-2u - p)K^{u - p)^ , 

where P^ 2 ^ is the completely antisymmetric project: Pi 2 ^ = | (id — Pu)- Using the crossing 
relations of the K-matrices (see ()4.9j) - (l4.12j) below) and after a tedious calculation (for details 
see Appendix A), we find that the quantum determinants of the most general K-matrices 
K^{u) given by (ITBl) and respectively are 
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Substituting dsn-dsn into dSISD, one has that 
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For generic rj, the fusion hierarchy does not truncate (c.f. the roots of unity case in 
the trigonometric limit OCg). Hence constitute an inhnite hierarchy, namely, 

j taking values 1, • •• The commutativity (EnD of the fused transfer matrices 

and the fusion relation (HID imply that the corresponding eigenvalue of the transfer matrix 
denoted by satisfies the following hierarchy relation 

A^^\u + rj — 2jrf) = A^^~^\u — 2{j — ^)? 7 ) A(m) — 5{u) A^^~^\u — rj — 2{j — l)ri), 

j = (3.16) 

Here we have used the convention A{u) = and = 1. Dividing both sides of 

()3.16|) by — 2{j — l)r]), we have 


A^^\u + r] - 2jr]) A^^ - r] - 2{j - l)r]) 

= T7— TD - it: -fVT + 


(3.17) 


Ab 2 )(m - 2(j - i)r7) ' ' ' Ab l)(u-2{j - ^)7]) 

We now consider the limit j —> cx). We make the fundamental assumption (Q (in 
particular, that the limit exists), which implies for the corresponding eigenvalues 


Q{u) = lim A^^\u — 2jri), (3.18) 

j —>+oo 

where we have used the same notation for the operator Q and its eigenvalues (c.f. Eq. (ESI)). 
It follows from (ITT7I) that 


Q M + T] Q{u-v) 
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Assuming the function Q{u) has the decomposition Q{u) = f{u)Q{u) with 
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where M is certain non-negative integer and {uj} are some parameters which will be specihed 
later (see and ()5.19j) below), then Eq. (I3.19j) becomes 


A(m) = Hi{u) 


Qiu + vi) 
Q{u) 


+ H^{u) 


Q{u - r]) 
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(3.21) 


Here Hi{u) = and H 2 {u) = 6{u) If is easy to see that the functions {Hi{u)\i = 

1,2} satisfy the relation 


Hi{u - rj)H 2 {u) = 5{u) 


(3.22) 
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with 6(u) given by ()3.15|1 . 

In summary, the eigenvalue A(m) of the fundamental transfer matrix t(u) (12.1 dll has the 
decomposition form (EnD, where the coefficient functions {ffi(u)} satisfy the constraint 
(Km . In the following, we shall use certain properties of the eigenvalue A(m) derived from 
the transfer matrix to determine the functions {Hi{u)} and therefore the eigenvalue A{u). 


4 Properties of the fundamental transfer matrix 


Here we will derive hve properties of the fundamental transfer matrix t{u), which together 
with the T-Q relation (II2H) enable us to determine the functions {Hi{u)}. 

In addition to the Riemann identity (jZl, the cr-function enjoys the following identity: 

2a{u)a{u + ^)a{u + ^)a{u — | f) 


a{2u) = 


1 T' 

’2 2 ' 


(j(m + 1 ) = —a{u), (j(m + r) = —V('u). 
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Moreover, from the definition of the elliptic function one may show that 
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M,2r), 

(4.3) 

M,2r). 

(4.4) 


The property iH) of the (T-function implies that the matrix elements of the K-matrices 
K^{u) given by (j2.15jl - (j2.16jl are analytic functions of u. Thus the commutativity of the 
transfer matrix t{u) and the analyticity of the R-matrix and K-matrices lead to the following 
analytic property: 


The eigenvalue of t{u) is an analytic function of u at finite u. (4-5) 

The quasi-periodicity of the elliptic functions, (|OI)-(|131), allows one to derive the following 
properties of the R-matrix and K-matrices: 


Ri,2(m + 1) = -0'i-Ri,2(M)aj'=-cr2Ri,2(M)cr2, i^u + 1) = - a"" ((o) a"", 

Ri2(M + r) = =-e"^“(“+2+i)a^R^,2(M)t^2. 

K-{u + t) = 

K+{u + t) = 
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From these relations one obtains the qnasi-periodic properties of t{u), 


t{u + 1) = t{u), t{u + r) = e-2*-(^+3)(2«+r?+r) ( 4 _g) 

The initial conditions of the R-matrix -Ri2(0) = P 12 and K-matrices: -R“(0) = id, trK~^{0) = 
a{2ri)/a{ri), imply that the initial condition of the fnndamental transfer matrix t{u) is given 
by 

^(“)L=o = X id. (4.7) 

The identity 
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snggests the semi-classical property of the R-matrix lim^^o {c’'(h).Ri 2 (w)} = (t(-u) id. This 
enables ns to derive the semi-classical property of t{u), 

\im.{a‘^^{ri)t{u)}=a‘^^{u)\im{tr{K'^{u)K~{u))} x id. (4.8) 

rj^O'' ^ rj^O ^ ^ 

Now, we derive the crossing relation of t{u). For this pnrpose, we introdnce 

R-(m) = tr2 (Pi2Pi2(-2n - 2p){iF2'(w)}*") , (4.9) 

Kt{u) = P 2 (Pi2Pi2(2n){P2^(n)}*=) . (4.10) 

Throngh a straightforward calcnlation (see Appendix A for details), we hnd 

K-{u) = -^^^ayK-{-u-r])ay = f4u)VK-{-u-r])V, (4.11) 

a{r]) 

K~^{u) = — — ^^Y^ayK~^{—u — ri)ay = f+{u)VK^{—u — ri)V, (4.12) 

where the 2x2 matrix V is eiven in (12.101). /_(«) = ^ and U(n) = The 


PT-symmetry (EZD and crossing relation (ITTUD of the R-matrix implies 
Tq°{-U-T]) = {Ro^Ni-u-T]) ...Ro^i{-u-T])y° 

= {Vo {RoAu )... RoMY^-'" Voy° 

= Ro {7?o,w(w)... Po,i(w)}*“*^-*" K) 

(-1)^-1 VoRo,i{u )... R^,n{u) Vo ^ (-l)"^"Vofo(n)Ro, (4.13) 
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Similarly, we have another “duality relations” between the monodromoy matrices T{u) and 

f(«), 

= (-l)''-‘V'„r„(u) V„. (4.14) 


The crossing relations of the K-matrices (j4.1 Ijl and (j4.1 and the “duality relations” (j4.1 djl 
and (j4.14|) . imply that 


t{—u — Tj) 


tro { {-u - r])To{-u - r])KQ {-u - r])To{-u 

tro 





\to 


-oi-u - v)Toi-u - ri)Y° jitTo {-u - ri)fo{-u - 
tro - v)iKoi-u - v)) 

tro (fo{u) {Vo{K+{-u - v)y°Vo} To{u) {Vo{K^{-u - v)Y°Vo}) 
tr„ (Ti(«)(A'+('t.))‘”T„(«)(A7(«))'") //+(«)/_(«) 
lOJ.llIinj trotritr2(^o{u)PoiRo,i{‘^u)K+{u)To{u)Po2Ro,2{-‘^u - 2ri)K:^{u 




f+{u)f-{u) 

trotritr2[Poifi{u)Ro,i{2u)To{u)Po2Ro,2{-2,u-2ri)K:^{u)Ky{u 


f+{u)f-{u) 

trotritr2[Ky{u)Ti{u)PoiRo,i{2u)Po2Ro,2{-‘2u-27])K2 {u)fi{u 


, 4 w 4 4 • (4.15) 

In deriving the second last equality of the above equation, we have used the relation, 
Ti{u)Ro,i{2u)To{u) = To{u)Ro^i{2u)Ti{u), which is a simple consequence of the quantum 
Yang-Baxter equation ()2.f)ll . Then, let us consider 


trotr2 {PoiRo,i{2u)Po2Ro,2(,-2u - 2r])K^ (u)) /f+(u)f_(u) 

= tro (PoiRo,i( 2 u)tr 2 {Po 2 Ro, 2 (- 2 u - 277)7^2"(«)}) /f+(u)f-(u) 

^ fro (^PoiRo,i{2u) //+(«)/-(«) 

^ tro (VoPoiRo,i{2n)Vo {K^{-u - 77 )}*°) //+(«) 

^ tro {ViPoiRoA2u)Vi {K^{-u - 77 )}*°) //+(«) 

= ViK^ {-u - r])Vi/f+{u) = ^ ^ [u) = {u) 

Substituting the above equation into (I4.15|l . we have the following crossing relation of t{u): 

t{—u — rj) = tri {u)Ti{u)K^{u)Ti{u'Yj = t{u). (4.16) 
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We remark that our above proof of the crossing relation (ITOll generalizes the previous proofs 
given in [ig. 


5 Eigenvalues and Bethe Ansatz equations 

It follows from (ji3I)-(I01) and (lorni that the eigenvalue A(n), as a function of u, has the 
following properties, 


Analyticity : A(m) is an analytic function of u at hnite u. 

Quasi-Periodicity : k{u + 1) = A(n), A(n + r) = e-2*-(^+3)(2«+^+PA(M), 


Crossing symmetry : 
Initial condition : 


A(-m 

A(0) 


-p) = A(n), 

a(p) ’ 


Semi-classic property : 


lira{tr{K~^{u)K («))} . 

rj ^ O '' ^ 77^0 ^ 


(5.1) 

(5.2) 

(5.3) 

(5.4) 

(5.5) 


TI 16 T-Q rGlcitioiis fj3.2ip ciiici dsm, together with the above properties dSIIll-dESl), can be 
used to determine {Hi{u)} and the eigenvalues of the fundamental transfer matrix. 

For this purpose, we dehne 6 discrete parameters = ±1|7 = ±, / = 1,2,3} and £x 
= -1-1 in the following. Associated with we introduce 


»F(«|{7’}) 




(T^^(u)cr( 2 u) TT A *7(7* + *? 

0-2^ (m r])a{2u + 2r]) TT iq 

a^^{r])a{2u + rj) f=T cT(e{^^a{^^) 


(5.6) 

(5.7) 


One may readily check that the functions H['^\u\{ep'^}) and H 2 '^\u\{ep'^}) indeed satisfy 
(KT2\i . namely. 


- ')l{P})J?f AKP}) = '5(«). (5.8) 

The general solution to can be written as follows: 

Hi(u) = i/!*>(«|{P})Si(«). ^^ 2 («) = fff>(«|{T}) 9 . 2 (f), (5.9) 

where {gi{u)} satisfy the following relations, 

9 i{u-v)g 2 iu) = 1, gi{u+1) = gi{u), g 2 iu + 1) = g 2 {u). (5.10) 
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(5.11) 


The solutions to (I^TTUl) have the following form, 

nfii - ut) 


1 nS. -U-- ,,) 


9i{u) = A: 


92{u) = 


nfi <T(u-ujy A nfi. <t(« -«+ -,,) ’ 

where A^i and N 2 are integers such that Ni, N 2 > 0, and A is an non-zero constant. In the 
above equation, we assume that uj 7 ^ — rj and u'j' 7 ^ '^ap ^ — 9 , otherwise the 

corresponding factors in gi{u) make transitions among 

respectively). Then the analyticity of A(-u) (I5.1|l requires that 9 i{u) and 92 {u) have common 
poles; i.e., iVi = iV 2 , and uJ = up + 9 . This means 

9 ) 


f \ X TT 
c/i(m) = AJ_J_ 

t=i 


a[u 


Ua 


92{u) = jI[ 


1 ^ o-[u 


Ua 


i=i 


a[u 


Ua) 


(5.12) 


Since HP^\u\{epp) = Hp^\—u — 9 \{epp), the crossing symmetry of A(n) (I5.d|l implies that 
H 2 {u) = Hi{—u — 9 ). Hence, A^i is even and 


2 


9i{u) = 


92{u) = 


i=i 


Ex 

2 




a{u - Uj + 9 )a{u + Uj + 29 ) 
a{u — uJ )a{u + uj + 9 ) 

9 )a{u + u~) 


a[u 


Ua 




A = ±1. 


(5.13) 


(5.14) 


This is equivalent to having additional Bethe roots; and the corresponding factors, except 
A, can be absorbed into those of Q{u) ()3.2()j) . Moreover the initial condition D implies 
A = -|-1. Therefore the eigenvalues of the transfer matrix take the following forms: 


(±)/ 1/^(7)! A T h) I Tj{±) j I r^( 7 )lA Qi'^ V) 


Hrpn\{eP’}) 


+ Hrpn\{eP>}) 


(5.15) 


Q{U) ^ V • 

In the above expression, we still have an unknown non-negative integer M defined in 
to be hxed. Moreover, A(m) must fulhll the properties (j5.2jl and (j5.5j) . This provides the 
following constraint to the boundary parameters 

3 3 


Y niod(l), JJ JJ 

7 =± l=l 


A')') _ 


= - 1 , 


(5.16) 


7 =± 1=1 

where k is an integer such that \k\ < N — 1 and N — 1 + k being even, and yields two different 
values for M in ()3.20|1 . corresponding respectively to Hj;'^\u\{epp), 


= -{N -iTk) 


(5.17) 
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The proof of the above constraint is relegated to Appendix B. 

Finally, we have that if the boundary parameters satisfy any of the constraints 

()5.16|) . the eigenvalues of the fundamental transfer matrix t{u) are given by 










- 7 ]) 


(5.18) 


where ^ ~ u^^'^)a{u + + rj) and the parameters respectively 

satisfy the Bethe Ansatz equations, 


ffhAf'uA}) 




(5.19) 


Indeed one can verify that both A^^'^{u) given by (I^TTHl) have the desirable properties (EH)- 
(PxHll provided that the constraint (j5.1b|l and the Bethe Ansatz equations (j5.19|l are satished. 
In the trigonometric limit r —> +icx), we recover the results in [H] (for details see Appendix 
C). The completeness [IZlini of the eigenvalues {A(^)(m)} in the trigonometric case suggests 
that for a given set of bulk and boundary parameters satisfying the constraint (I5.1b|l . the 
eigenvalues A*^“)(m) and together constitute the complete set of eigenvalues of the 

transfer matrix t{u) of the open XYZ spin chain. 

Therefore the complete set of the energy eigenvalues of the Hamiltonian (jl.l|l . when 
the boundary parameters satisfy the constraint (I5.16|) . are respectively given by 


( \ ^ ^ 

2 ^ (c(uf>) - C(«f + ’;))+(^ -+ E . (5.20) 

^ ^ [ J=1 7=± «=1 J 

where the parameters respectively satisfy the Bethe Ansatz equations (Em- 

Some remarks are in order. Firstly, in |TH] a generalized algebraic Bethe Ansatz [TH] 
was used to diagonalize the open XYZ spin chain when boundary parameters satisfying 
a constraint. However, only partial eigenvalues of transfer matrix corresponding to our 
A^~^\u) were obtained there and it is not clear whether the approach in ^H] can give all the 
eigenvalues of the transfer matrix. (See also (20]; where partial eigenvalues of the transfer 
matrix of the open XXZ spin chain were obtained.) Indeed, it is found in this paper that 
both A(^(’(m) are needed to constitute a complete set of eigenvalues. This implies that 
two sets of Bethe Ansatz equation (ICTl) and therefore presumably two pseudo-vacua are 
required. Unfortunately, it is not yet clear how to construct the second pseudo-vacuum in 
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the generalized algebraic Bethe Ansatz approach. Secondly, when the crossing parameter rj is 
equal to with p being a non-negative integer (i.e. roots of unity case), the corresponding 
fusion hierarchy truncates as in the trigonometric case ^3], and one may generalize the 
method developed in ra to obtain the corresponding functional relations obeyed by the 
transfer matrices and therefore the eigenvalues of the fundamental transfer matrix t{u). If 
the boundary parameters satisfy the constraint (15.161) . the complete set of eigenvalues and 
the associated Bethe Ansatz equations, in the roots of unity case, are expected to be still 
given by and respectively. Thirdly, it would be interesting to determine the 

conditions for which the limit dn exists. We have seen that (IK:T^ - (ICT1 solve the open 
XYZ chain with generic values of rj if the constraint (I5.16|) is satished. This suggests that, 
for generic values of p, the constraint ()5.16|) may be a necessary condition for the existence 
of the limit dH. It would be interesting to explicitly evaluate the Q-operator directly from 
Eq. dH- 

6 Conclusions 

We have argued that the Baxter’s Q-operator for the open XYZ chain is given by the j —> cxo 
limit of the transfer matrix with spin-j auxiliary space. This together with the fusion hier¬ 
archy leads to the Baxter’s T-Q relations and (IT^ . These T-Q relations, together 

with the additional properties dsm-dEsi), enable us to successfully determine the eigenval¬ 
ues flOHD of the fundamental transfer matrix and the associated Bethe Ansatz equations 
onni). For a given set of bulk and boundary parameters satisfying the constraint ()5.16|) . 
the eigenvalues and together constitute the complete set of eigenvalues of 

the fundamental transfer matrix t{u), which leads to the complete spectrum (I5.2()j) of the 
Hamiltonian dm. Our derivation for the Q-operator and associated T-Q relation can be 
applied to other models which share s/ 2 -hke fusion rule. Moreover, it would be interesting to 
generalize our method to integrable models associated with higher rank algebras Eni^am. 
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Appendix A: Determinants and crossing relations 


In this appendix, we compnte the determinants (ITTTl) and (ITT^ of the K-matrices. In so 
doing, we will also prove the crossing relations Km and (I4l2l) . 

One rewrites K~{u) in (HH as 

3 

K~{u) = ^ fa{u)a°‘, cr° = I, = a^, = a^, 

a=0 

It is easy checked that the coefficient fnnction {/o(m)} satisfy 


foi-u) = /o(m), fi{-u) = -fi{u), for i = 1, 2, 3. 


The above parity relations lead to the following nnitarity relation of K (u) 

3 

K~{u)K~{-u) = A 2 {u) X id, A 2 (n) = ^ 

0=0 


From the expression of K (u) and after direct calculation, we hnd that 


A2(m) = f[ ^ ^ ^ 






1=1 


(A.l) 


(A. 2 ) 


We now follow a method similar to that developed in [21]. Noting -Ri, 2 (—h) = ^ 

PT-symmetry (EH) of the R-matrix, the reflection equation (ITTTl) . when ui = u, U 2 = u + T], 
reads 


Fi 2 (m)Fi^ 2 ( 2 m + ifjK^ {u + ini) = {u + rf)Ri^ 2 {‘^u + r])K^ {u)Pi 2 
Since the dimension of the image of P 12 ^ is equal to one, the above equation becomes 
K2{u + r])Ri^2{‘2u + r])K^{u)Pi2'‘ = +'^)-R2'('“ + h)-Pl2'^ 

= tri 2 (^Pu ^Ki (m)Fi, 2 ( 2 m + r])K 2 {u + 77 ) j 

(A.3) 


Det{K-{u + r])}Pi2^. 


It is well-known that the completely anti-symmetric vector with unit normal in ® is 

2 


^ (|1) ® |2) - |2) 0 |1)) = (id 0 V) N) ® N) ) > 


(A.4) 
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where {|1), |2)} are the orthnormal basis of C^. Acting both sides of on the vector 

()A.4fl and after straightforward calcnlating, we hnd that is indeed eqnivalent to 

K-{u) = Det{K-{u + r])} {VK-{u + . 

Using (EH), we have 


^ Det{A: {u + T])} 

K (u) = -T—;-r- VK [—u — riW. 

^ ’ A 2 (n + r 7 ) ^ 

Taking trace of both sides of ()A.5jl . we have, for the L.H.S. of the resnlting relation, 

L.H.S. = tr {K-{u)) ® tri 2 {PuRi,2{-‘2u - 2r]){K-{u)y^) 

= tr2 {tri {Pi2-Ri,2(-2n - 27])} {K^ (m)}*^) 

= {a{—2u — 2r]) + c{—2u — 2//)) tr 2 [K^iu)) 

a(2u) 


(A.5) 


= {a{—2u — 2rj) + c{—2u — 2rj)) 
(j{2u)(j{2u + 27]) 


a{u) 


(A.6) 


a{r])a{u + t]) 

In the above deriving, we have used the following relations: 

tri (Pi 2 Pi, 2 (m)) = (aiu) + c{u)) x id, 

/ N / N a(2u + 2r])a{u) 

a{-2u - 2r]) + c{-2u - 27]) = - ( \ ( . \ ^ 

a{7])a{u + 7]) 

which can be obtained from (Q-(E3D, and the following identities of the elliptic functions: 


- 1 - 
2 

1 

(2m, 2t) 

= e 

- 1 - 
2 

1 

- 2 - 



- 2 - 

■ 0 ■ 
1 

- 2 . 

(2m, 2r) 

= e 

■ 0 ■ 
1 

. 2 . 


(r.2r) X + 5 


(0,2r) X 


a{u - ^)a{u + I + |) 


On the other hand, for the R.H.S. of the resulting relation, we obtain 

Det{A'“(n + t])} 


R.H.S. = 


A2{u + 7]) 
Det{K~{u + 7])} 
A2{u + 7]) 
Det{K~{u + 7])} 


tr {VK-{-u- 7])V) 
tr (y‘^K~(—u — 7])) 
tr (^K~{—u — T])') 


A2{u + 7]) 

Det{K~{u + 7])} a{2u + 2p) 


A2 {u + 7]) a{u + 7]) 


(A.7) 
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Comparing with 


3)) and (EH), we have 

Det{iC“(M + 77 )} a{2u) 


(A.8) 


A2{u + ri) a{r]) ' 

Now (jd.ld|) follows from ()A.2j) and ()A.8|) . Similarly one can prove (I8.14j) using (I2.1fi|) . Finally, 
(IAaIi and ()A.8|1 give rise to dHH), and (ICT) can be similarly be proven from (ETTBD . 


Appendix B: Constraint condition 


One can easily check that both forms of eigenvalues of t{u) satisfy the hrst quasi¬ 

periodicity of (EH as required. Now we investigate the second quasi-periodic property. The 
dehnitions (j5.b|l and ()5.7|1 of imply 

Keeping the above quasi-periodic properties and the dehnition ()8.2()j) of Q{u) in mind, in 
order that have the second property of (Q one needs that the boundary parameters 

satisfy the following constraint 

3 

= kr] mod(l), (B.l) 

7 =± l=l 

where k is an integer such that \k\ < N — 1 and N — 1 + k being even. For such integer k, let 
us introduce two non-negative integers | (N — 1 =p /c). Moreover, associated 

respectively with one introduces Q^^\u), 

m(±) 

Q^^\u) = a{u - uf"^)a{u + uf"'^ + rj). 
j=i 

Hence the eigenvalues of the transfer matrix take the following forms: 

+ AAukA}) (B.2) 

At this stage, we still have freedom to chose the discrete parameters Now we apply 

the semi-classical property (E31) of A(m) to obtain a restriction to this freedom. 

For this purpose, let us introduce A(^)(m), corresponding to (El, 


A(±)(u) = 




gw(M) ’ 


(B.3) 
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and define = lim^^o ■ Since the value of the discrete parameters does not 

change when taking the limit of 77 ^ 0 , the constraint (EH) in this limit becomes 

° mod(l). (B.4) 

7=± 1=1 

Substituting (El into the semi-classical property dESD of the eigenvalues of the fundamental 
transfer matrix t{u), we have for the L.H.S. of the resulting relation, 


L.H.S. = lirn {cr^^(? 7 )A^^^(M)} 


= 


nn 




,(7);:r,(7)^ 


. 7 =± 1=1 


TT TT q'(-» =F 

1111 

7=± i=l ) 


Comparing both side of the resulting relation, (jEH is equivalent to the following relation 


3 / , ( 7 )--( 7 )\ 

M a{u ± 6 ; ai ) 

7=± i=l ) 


lim{tKA-WA-W)} . (B.5) 


7=± 1=1 ) 

Keeping the constraint (jR.4j) among after straightforward calculation, we hnd that 

(El is satished if and only if the following constraint among the discrete parameters 
were fulhlled 

3 


nn 


,( 7 ) 


= -1. 


(B. 6 ) 


7 =± 1=1 

(El and ()B.6|1 are nothing but the constraints listed in (15 .161) . 


Appendix C: Trigonometric limit 

In this Appendix, we consider the trigonometric limit r —>■ +ioo of our result. The dehnition 
of the elliptic functions (EU-dZl imply 


a{u+^) = 9 


0 
1 

L 2 J 




(C.l) 


and the following asymptotic behaviors 


lim a{u) = —2e 4 sin nu -|- 

r^+zoo 


lim 9 

r^+zcxD 


0 

1 

L 2 


{u,t) = 1 + .... 
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(C.2) 

(C.3) 











The above asymptotic behaviors lead to the well-known result: 


lim R{u) 

r^+ioo 


1 


/ sin 7i{u + rj) 


sin TTrj 


\ 


\ 

sin TTU sin nrj 
sin Tir] sin ttu 

sin niu + t]) J 


After reparameterizing the boundary parameters 






a 


(T) _ 


= TA 


^ + 2’ 


OL 


(T) 


= -B 


i T 
^ 2 2 ’ 


and using (Unii-ra, one has 


lim K (u) = Kf,-. (u) = -— |2 sin vra. cos 7i(5_ cos nu I 

T^+ioo '■o 2 sin 7rQ;_ cos 7rp_ 

-|- 2i cos TTtt- sin tt(5- sin ttu -|- sin 2?™ (cos ct* — sin a^)} , 


lini a:+(m) = a:J)(“) = 


^+200 


^{t)\ 


-u 


v) 


ia-,0-,9-)^{-a+,-P+,e+) 


(C.4) 


(C.5) 


(C.6) 

(C.7) 


Hence our R-matrix and K-matrices, when taking the trigonometric limit, reduce to the cor¬ 
responding trigonometric ones used in j3] after a rescaling of the corresponding parameters. 
In order for the constraint ()5.16j) to be still satished after the limit, one needs further to 
require that = —63"'"^ Then, the resulting constraint becomes that in [^. Moreover, our 
in the limit give rise to the corresponding trigonometric ones in up to 
some constants (due to the constant factors appearing in ()(k4ll . (1(1.611 - (1(1.71) 1. Therefore, our 
result recovers that of [S] in the trigonometric limit. 
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